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The cross-spectral acoustic analogy is used to predict auto-spectra and cross-spectra 
of broadband shock-associated noise in the near-field and far-field from a range of heated 
and unheated supersonic off-design jets. A single equivalent source model is proposed for 
the near-field, mid-field, and far-field terms, that contains flow-field statistics of the shock 
wave shear layer interactions. Flow-field statistics are modeled based upon experimental 
observation and computational fluid dynamics solutions. An axisymmetric assumption is 
used to reduce the model to a closed-form equation involving a double summation over 
the equivalent source at each shock wave shear layer interaction. Predictions are compared 
with a wide variety of measurements at numerous jet Mach numbers and temperature ratios 
from multiple facilities. Auto-spectral predictions of broadband shock-associated noise in 
the near-field and far-field capture trends observed in measurement and other prediction 
theories. Predictions of spatial coherence of broadband shock-associated noise accurately 
capture the peak coherent intensity, frequency, and spectral width. 


Nomenclature 


Symbols 

Description 

Uj 

Fully-expanded jet velocity 

A ijt m 

Coefficient matrix 

X 

Observer vector 

D 

Nozzle exit diameter 

y 

Source vector 

Dj 

Fully-expanded diameter 

Vc 

Jet core length 

c 

Speed of sound 

z 

Radial cylindrical coordinate 

Ft 

Far-field term 

Oti 

Constant 

fi 

Equivalent source term 

p 

Off-design parameter 

G 

Cross-power spectral density 

r 

Coherence 

I 

Modified Bessel function 

7 

Ratio of specific heats 

k 

Turbulent kinetic energy 

A 

Separation operator 

la 

Turbulent length scale coefficient 

v(€>ri,C) 

Source separation vector 

u 

Turbulent length scale 

P 

Density 

M d 

Design Mach number 

Qi 

Observer angle relative 

M t 

Mid-field term 


to downstream nozzle axis 

Mj 

Fully-expanded Mach number 

<Pi 

Azimuthal cylindrical coordinate 

N s 

Number of shocks 

V 

Observer angle relative 

N t 

Near-field term 


to upstream nozzle axis 

Pf 

Constant 

T 

Retarded time 

P 

Pressure 

T c 

Turbulent time scale coefficient 

Ps 

Shock pressure 

T s 

Turbulent time scale 

Pt 

Total pressure 

UJ 

Radial frequency 

R 

Gas constant or radiation distance 

UJ 

Modified radial frequency 

R mn 

Two-point cross-correlation 




Two-point cross-correlation of source 

Abbreviations 


Si 

Shock wave thickness 

BBSAN 

Broadband shock-associated noise 

St 

Strouhal number 

CPSD 

Cross-power spectral density 

r 

Radiation vector 

PSD 

Power spectral density 

T 

Temperature 

RANS 

Reynolds-averaged Navier-Stokes 

Tj 

Fully-expanded temperature 

SHJAR 

Small Hot Jet Acoustic Rig 

Tij 

Lighthill stress tensor 

TTR 

Total temperature ratio 

u 

Velocity vector 



u c 

Convection velocity 




*Research Aerospace Engineer, Aeroacoustics Branch, NASA Langley Research Center, 2 N. Dryden St. MS 461, Hampton, 
VA, 23681, USA, AIAA Senior Member, s.miller@nasa.gov 


1 of 21 


American Institute of Aeronautics and Astronautics 


Introduction 


The interaction of high speed convecting turbulence with a shock wave results in high intensity acoustic 
radiation. This radiation is broadband in nature, and it is dependent on the convection velocity and intensity 
of the turbulence, and strength of the shock wave. Systems of shock waves are present within off-design 
supersonic turbulent jets. Off-design supersonic jets are created by propulsion systems of flight vehicles 
that include rockets, missiles, aircraft, and by some natural phenomenon. Coherent turbulence in the jet 
shear layer interacts with the shock cells and produces broadband shock-associated noise (BBSAN). The 
radiation from subsequent shock wave shear layer interactions combines constructively and is characterized 
by multiple broad spectral lobes. This radiation may impinge on the flight vehicle and can also cause observer 
hearing loss or annoyance. Predicting the cross-power spectral density (CPSD) of BBSAN is beneficial for 
understanding the physics of the radiating source and also provides mathematical guidance for flight vehicle 
designers. This paper presents an approach to predict the near-field intensity and coherence of off-design 
supersonic jets using the cross-spectral acoustic analogy. 

In the far-field BBSAN is dominant at mid- to high- frequencies in the upstream and sideline direction 
relative to the jet flow. Generally, BBSAN is dominated by other noise components in the downstream 
direction. Harper-Bourne and Fisher 1 proposed that the intensity of BBSAN increases as /3 4 , where /3 is 
the off-design parameter involving the fully-expanded Mach number, Mj, and design Mach number, M ( j : . 
A general form of the off-design parameter is, /? = M'j — Mjl 1 / 2 . Norum and Seiner 2 showed that the /3 4 
relationship is invalid when Mach disks form within the jet plume with increasing nozzle pressure ratio (NPR) . 
Generally, the intensity scaling of BBSAN is strongly dependent on NPR and very weakly dependent on total 
temperature ratio (TTR). Viswanathan et al. 3 noted that when NPR is held constant and TTR increases 
the BBSAN intensity saturates (ceases to increase). A physical explanation for the saturation of BBSAN was 
proposed by Miller. 4 The source of BBSAN is the interaction of large-scale coherent turbulence within the 
jet shear layer interacting with the shock waves. Each of these shock wave shear layer interactions creates 
outgoing acoustic waves. As the waves propagate from the near-field to the far-field, they constructively 
interfere to form characteristic broad lobes in the far-field. BBSAN spectra are characterized by multiple 
broad lobes that decrease in intensity with increasing frequency. Though this noise component is broadband 
and random, it is highly coherent near peak frequencies due to its constructive nature, and almost entirely 
incoherent at relatively low and high frequencies. The spatial coherence of BBSAN decreases with increasing 
observer separation distance and increasing frequency. 

A number of experiments were conducted to study the noise from off-design jets. Early measurements of 
Yu 5 yielded contour plots of BBSAN intensity in the jet near-field. Soon after, Tanna 6 measured far-field 
auto-spectra for a wide range of jet Mach numbers and temperature ratios from convergent nozzles, and 
isolated the BBSAN component of the spectrum using the prediction method of Harper-Bourne and Fisher. 1 
Tam and Tanna' defined /3 to account for convergent-divergent nozzles and jet stagnation temperature 
variation. The same year that Tam and Tanna' presented their work, Norum and Seiner 2 examined the 
variation of peak Helmholtz number of BBSAN relative to unheated jets from a convergent nozzle. Norum 
and Seiner 2 used a near-field microphone array to characterize BBSAN intensity and their results agreed with 
the early findings of Yu. 5 Seiner and Yu 8 used a near-field linear microphone array located 2.68 jet diameters 
from the nozzle centerline to confirm measurements of Harper-Bourne and Fisher, 1 who showed that the 
shock-associated noise intensity is proportional to the first power of the fluctuating velocity components 
within the jet mixing layer. 

Bridges and Brown 9 and Bridges and Wernet 10 examined far-field auto-spectra and performed time 
resolved particle image velocimetry from a large range of off-design jets. Turbulent statistics along the nozzle 
lip line were shown to possess some similarity for on-design jets. For off-design jets, the similarity is not 
readily apparent because of the large gradients of turbulent kinetic energy (k) (and other turbulence statistics) 
after each shock wave shear layer interaction. The far-field auto-spectral and turbulence measurements of 
Bridges and Brown 9 and Bridges and Wernet 10 are an important part of the present model development. 
These increases in turbulent kinetic energy are complicated by oscillations of the shock waves (as described 
by Panda 11 ) and cause a noticeable effect of broadening the spectrum of BBSAN in the far-field. 

Recently, Viswanathan et al. 3 isolated the BBSAN component of jet noise incoherently, like Tanna, 6 but 
used the theory of fine- and large-scale similarity spectral subtraction. They observed nonlinear propagation 
of BBSAN, the saturation of BBSAN with increasing temperature, and most importantly the decay of 
coherence between observer pairs. A similar investigation of Kuo et al., 12 using helium-air mixtures to 
simulate jet heating, showed that BBSAN saturates with increasing temperature. Savarese et al. 13 made 
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measurements within the jet near-field for multiple pressure ratios. Significant coherence was observed 
between fluctuating velocity quantities at the shock wave shear layer interactions and the acoustic field. 
Savarese et al . 13 also examined flight vehicle Mach number effects and their measurements were in agreement 
with Norum and Shearin 14 and Ahuja et al . 15 

Many mathematical models were developed to predict BBSAN, and here we survey those without signifi- 
cant empiricism. Harper-Bourne and Fisher 1 created the first successful prediction methodology for BBSAN 
that depends on the rate of decay of turbulence between shocks and a characteristic spectrum produced 
by each shock wave shear layer interaction. The model was based on the premise of coherent interaction 
between turbulence in the jet shear layer and the nearly periodic jet shock cell structure. More recently, 
Harper-Bourne 16 created a prediction technique for the near-field intensity of jet mixing noise and BB- 
SAN. This near-field model for BBSAN intensity depends on a master spectrum, which is like the model of 
Harper-Bourne and Fisher . 1 

Tam et al., 1 ' based on the preliminary work of Tam and Jackson , 18 developed a linear multiple scales 
model of the shock-cell structure. This model, which is more physical than the model of Pack , 19 takes into 
account the gradual spatial change of the jet spreading and the smoothing of gradients by turbulence. The 
effect of divergence and turbulent dissipation was included following the method of Tam et al. 1 ' Tam 20 
developed a stochastic model for BBSAN where the basic physical model was described by Tam and Tanna. ' 
The large-scale turbulence in the jet shear layer was modeled as a random superposition of instability waves 
supported by the jet mean flow, as described by Tam and Chen . 21 Tam 22 modified the model of Tam 20 to 
predict BBSAN from heated jets with a moderate off-design parameter. The stochastic model of Tam 20 was 
shown to predict near-field BBSAN intensity and compared very favorably with the measurements of Yu . 5 

Morris and Miller 23 developed a prediction method for BBSAN that is based upon an acoustic analogy 
approach. The Euler equations were rearranged into the linearized Euler equations (LEE) operator and 
were equated to equivalent sources. The resultant model consisted of a volume integral containing the jet 
plume and an integral of the shock cell pressure wavenumber spectrum. Their acoustic analogy successfully 
predicted BBSAN for a wide range of jet operating conditions. Miller and Morris 24 studied propagation 
of BBSAN by altering the vector Green’s function of the LEE to contain the Green’s function of Lilley’s 
equation as an argument. It was shown that in the upstream and sideline radiation direction, refraction 
effects have little impact on BBSAN. Very recently, Miller 25 created an acoustic analogy for jet noise that 
included equivalent sources for both turbulent mixing noise and BBSAN, that was shown to predict the total 
noise accurately for a wide range of jet Mach numbers and jet stagnation temperatures. 

The models of Harper-Bourne and Fisher , 1 Harper-Bourne , 16 Tam , 20 or Miller 25 are very good at pre- 
dicting the power spectral density in the far-field; however, they are not capable of predicting cross-power 
spectral density in the near-field. Recently, the cross-spectral acoustic analogy was developed by Miller 26 
to overcome similar difficulties present within other problems of interest to NASA. It can be shown that 
the cross-spectral acoustic analogy can be simplified to produce Lighthill’s 27, 28 acoustic analogy when both 
observers are in the far-field and at the same location. 

In this paper we formulate a model for the CPSD of BBSAN based upon the cross-spectral acoustic 
analogy of Miller . 26 A source term is proposed for the two-point cross-correlation of the Lighthill stress 
tensor and is valid for a wide range of jet operating conditions. Arguments of the source term are dependent 
on the jet Mach number and jet temperature, among other terms. The scaling of the source model is 
formulated based upon the work of Morris and Miller 23 but in the framework of Lighthill’s acoustic analogy. 
Anisotropic effects of the source are partially taken into account through arguments of the observer vectors 
and varying turbulent scales. Spatial integration of the shock wave shear layer interactions are performed 
analytically by assuming that the jet flow-field is axisymmetric. The resulting closed-form model for the 
CPSD of BBSAN in the near-field involves a double summation of the source at all shock wave shear layer 
interactions. 

The next section of this paper presents a mathematical model for the prediction of the CPSD of BBSAN. 
Arguments of the CPSD BBSAN model are defined in relation to the jet Mach number and temperature. 
Predictions are presented for far-field and near-field BBSAN intensity and BBSAN spatial coherence and are 
compared with measurement. Finally, a summary and conclusion closes the paper. 
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Mathematical Theory 


We assume that the Navier-Stokes equations model the physics of shock wave shear layer interactions 
and subsequent sound production and propagation. A partial derivative operator with respect to time and 
a divergence operator are applied to the continuity and momentum equations respectively. The difference 
of the resulting equations are added to a wave equation operator with mixed independent values of pressure 
and density. The resultant equation contains the well known Lighthill stress tensor, (see Lighthill 2 ' for 
details and discussion of T t j ) . The right hand side of this equation represents an equivalent aerodynamic 
source of sound. We expand the first and second terms of the right hand side, and unlike Lighthill, 2, we 
retain all terms. A vector free stream Mach number is introduced to account for source convection, density 
is converted to pressure, and the result is written more compactly. The resulting equation for pressure 
is used to form a cross-correlation between two observer points X\ and X 2 . The theories of Lighthill, 29,27 
Ribner, 30 and Ffowcs Williams 31 are used to simplify the resulting spatial cross-correlation for pressure using 
the assumption that the turbulence is statistically stationary. We then group integrand terms according to 
their relative contribution to the near-field, mid-field, and far-held. Terms with second order powers in 
propagation correspond to far-held terms, fourth order powers correspond to the mid-held terms, and sixth 
order powers correspond to near-held terms. The CPSD is then found using a forward transform. We hnd 
the cross-spectral acoustic analogy of Miller 26 


OO OO 


r Far-Field Term Mid-Field Term Near-Field Term ' 


G(x i , x 2 , uj ) = 


167T 2 


FtTijTl m + i m + NfTijT j,' 


•rim 


— OO — OO 


Ci) 


x exp 


■ , , n r 2 

—luj ( r -| 

Coo Cqo 


drdy 2 dy 1 , 


where c is the speed of sound, G is the cross-power spectral density, X\ and x 2 are vectors from the origin to 
observers, y 1 and y 2 are vectors from the origin to source positions, r is the retarded time, and w is the radial 
frequency. Vectors rq and r 2 are dehned from each source position to each observer. The primes denote 
quantities at position y 2 . Evaluation of Eqn. 1 yields the CPSD due to a convecting turbulent held within 
a free-stream Mach number, and its application is not limited to jet hows. Note that the propagation 

operator of the cross-spectral acoustic analogy is the three-dimensional wave equation operator. Two-point 
cross-correlations involving various retarded time derivatives of the Lighthill stress tensor, Tjj , appear before 
various terms. Prefactor terms, F tl M t , and N t , involve directional trigonometric functions that correspond 
to the far-held, mid-held, and near-held. For example 


Ft = 


' I' 3' l 1 m 


L C to rr ' 


( 2 ) 


and M t and N t are shown by Miller, 26 and are not repeated here for compactness. To evaluate Eqn. 1 
knowledge of TijT[ m and its derivatives must be known. In particular, for out-going acoustic radiation we 
must model TijT[ m and note the relation 


rp /in/ D V 

± % 3 ± lm — q t A ^ijlm 




( 3 ) 


where r/ = is a vector between source positions y x and y 2 . 

We must now model R)’ llrn (y-f , rj, r), and the model should characterize the BBSAN equivalent source. 
As there is no basis for an equivalent source for BBSAN in the context of a Lighthill acoustic analogy, we 
pause and reexamine the model of Morris and Miller 23 who developed a similar relation. The model of Morris 
and Miller 23 is developed for an acoustic analogy based on the linearized Euler equations. Their two-point 
cross-correlation of the right hand side source terms that characterize BBSAN is dehned as 


R V mn = I+V^ + t), 

where 


fi = -U S 


du t 


dxj 


u td 


du St j 
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and u s and ut involve unsteady velocity components due to the shocks and turbulence respectively. The 
term fV is present on the right hand side of the acoustic analogy shown in Morris and Miller. 23 Morris and 
Miller 23 proposed a model for the two-point cross-correlation of /)' as 


R v = 


PsPs, v k 
rp- rp /2 


exp 


exp 


(£ - U C UT ) 2 


l 2 



r (v 2 +c 2 )] 

exp 

[ 11 \ 


( 6 ) 


where k is the turbulent kinetic energy, p s is the shock pressure, l s and l± are streamwise and cross-stream 
turbulent length scales, u is the time-averaged streamwise velocity component, and t s is the turbulent time 
scale. Note that the subscript y denotes that the quantity is evaluated at y 1 + 17. Scales of turbulence vary 
spatially within the turbulent field, and can be approximated numerically with a steady Reynolds- averaged 
Navier-Stokes (RANS) solution, unsteady numerical methods, or with simpler models (that will be discussed 
below). Vector components of r) are £, 17, and £ respectively. 

Based on the model of Morris and Miller, 23 which was developed for the LEE operator and its associated 
equivalent sources, we propose an equivalent source for BBSAN 


where 
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Aijim is a coefficient matrix that describes anisotropic effects, and subscripts y and z denote cross-stream 
scales in radial directions. The terms chosen within Eqn. 8 are Guassian and model two-point cross- 
correlation measurements of velocity fluctuations within turbulent flows. Relatively more complicated forms 
or empirical fits can be chosen, but these choices do not yield large increases in prediction accuracy and pre- 
vent the ease of analytical integration. We scale the two-point cross-correlation with the local k. The terms 
of Eqn. 7 are chosen based on the model for turbulent mixing noise (see Miller 26 ) and Eqn. 4 . BBSAN inten- 
sity scales with NPR and TTR as approximately p 2 k , and for normalization and TTR saturation purposes 
is multiplied by c -2 . The scaling of BBSAN intensity using this equivalent source model compares favorably 
with / 3 4 . In summary this choice of R"ji m is consistent with the model of Morris and Miller 23 and Harper- 
Bourne and Fisher, 1 in that the scaling is consistent with their predictions and associated measurements, 
respectively. 

By examining the far- field, mid- field, and near- field terms of Eqn. 1 , it is apparent that Eqn. 7 must be 
differentiated with respect to r multiple times, multiplied by an exponential factor involving r, and finally 
integrated with respect with r. These operations for the near-field term result in 
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Note that the mid-field operations result in a negative sign. It is removed because the original CPSD 
formulation removed a corresponding mid-field sign. Thus, the formulations are consistent. Finally, the 
far-held operations result in 


& ( Aijim PsPs,r]k 


dr 4 1 7T 1 / 2 TgW 4 C 2 

Aijlm PsPs,r)k 
TT^Tfw 4 C 2 


exp 


H£l] 

exp 

L UT S J 



— (£ ~ UT ) 2 " 1 


P. 


exp 


_ 2 -| 


7 2 

sy 


exp 


A 

/ 2 


21 


exp 


-|£l 


4 ? ,4 

(3Z 4 - 12Z 2 (£ - St) 2 + 4(£ - Sr) 4 ) exp 


x exp 


l 2 

L sy J 


exp 


-c 2 


i 2 

l sz J 


exp [— iwr] dr 


-(g-nr) 2 

ll 


exp [—iuir] dr 


Aijlm PsPs,r}^ Is 


exp 


Hen 


r-r? 2 ] 


r-c 2 i 



exp 

k 2 J 

exp 

l p„ \ 

exp 

[ UT S J 





ui (4 iu^ + Z 2 w) 
4m 2 


( 11 ) 


The results of the operations shown in Eqns. 9 through 11 are now used in conjunction with Eqn. 1, and 
after simplifying 
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Equation 12 can be integrated numerically, but here we seek a simplified approach. We convert to 
cylindrical coordinates using the transform y 1 — > y 1 (yi, z\ cos[0i], z\ sin[0i]), where Z\ is the distance from 
axis y\ to the source. Similarly we use the relation y 2 — > y 2 {y 2 , z 2 cos[0 2 ], z 2 sin[0 2 ]). We assume that the 
how-held is axisymmetric about the centerline axis. Recall that the source of BBSAN is due to shock wave 
shear layer interactions. These interactions occur within very small volumes of space relative to the total 
volume of the jet plume. We assume that the source quantities vary very slowly spatially within the region 
of the shock wave shear layer interaction. These assumptions allow us to model the shock wave shear layer 
interactions as ‘rings’ of equivalent BBSAN sources centered about the nozzle axis. We note the geometric 
relation 
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Using these assumptions and Eqn. 13, Eqn. 12 is 
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(14) 


where Z sr . is the radial turbulent length scale, N s is the number of shock wave shear layer interactions, and S 
represents the local shock thickness in the axial (x) and radial (r) directions respectively. Like the coherence 
and CPSD of mixing noise, the mid-field and near-field terms will have characteristic peak frequencies lower 
than that due to the far-field term. Turbulent mixing noise intensity is dominant in the entire flow-field at 
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low frequencies. Thus, we focus our attention on the far-held term when examining the BBSAN component. 
This does not imply that the model is invalid in the near-held or mid-held because a far-held assumption 
has not been made. The integration involving (f> 2 and <p 1 within Eqn. 14 can now be conducted analytically. 
Because the shock wave shear layer interactions occur at relatively discrete locations compared to the total 
axial extent of the jet plume, the hnal axial integrals of £ and £' are approximated as a double summation 
involving the number of shock wave shear layer interactions. The hnal model equation for the CPSD of 
BBSAN is 
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where I o is the modihed Bessel function of the hrst kind of order zero. Note that Aiji m implies summation. 
Subscripts m and n imply evaluation at the rn th and n th shock location. Equation 15 is in closed-form 
with special function I 0 and can be used to predict the cross-spectrum from the shock wave shear layer 
interactions of axisymmetric jet hows. 

We must model the arguments of Eqn. 15 that are associated with the turbulence within the jet. The 
remaining terms need to be modeled. They could be found through numerical computation, but here we rely 
on empirical models generated from measurements. The core length of the jet, y c , is approximated using the 
equation of Eggers 32 and modihed by Tam et al. 17 


Vc_ = I 4.2 + 1.1 M] + 1.1(1 - Tj /Too) for Tj < 

D > 1 4.2 + 1.1 M] + exp[— 3.2(Tj/Too - 1)] - 1 for T. } > 

where Tj is the fully-expanded temperature. The number of shocks within the jet plume, N s , is estimated 
by performing the maximization 


N s := arg max 

N a 


Dc > Dj 



N s /3 + — — -N s 

sr 10 



(17) 


Equation 17 is based partly upon the theory of Pack 19 and calibrated by examining the number of shock 
wave shear layer interactions within a database of off-design jet steady RANS solutions. The fully-expanded 
diameter, D 3 . is adopted from Tam and Tanna. ' Equation 15 assumes sources are at relatively discrete 
locations. These locations are estimated using the theory of Pack, 19 and a simplifying assumption that they 
reside on the sonic line that extends roughly from the nozzle lip, = y^O, Dj/2, <j>), to approximately 
the end of the potential core, y i ~ y(y c ,0, 0). The shock pressure, y S) j, at each shock wave shear layer 
interaction, i = 2 to N s is based upon a linear estimation 


y Sii = -Ay + Ay + poo, (18) 

Vn„ 

where y Si i = Ay + for i = 1 and the isentropic formula 


Ay = p t 



Poo • 


(19) 


Here pt is the total pressure in the nozzle plenum. The streamwise velocity component is estimated using 
the empirical formula of Witze 33 


u = 
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for yi > y c , 


(20) 
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and convection velocities are estimated by multiplying u by 0.70. The variation of temperature, T, is assumed 
to possess the same form as the empirical relation for u. The variation of the speed of sound is c = ( yi?T ) 
and the variation of density is p = (p^ + p S j){RT)~ l . Scales of turbulence are estimated by 


(1 \ 

1 

r u> 2 11 

n D i + otiyiDj exp 

T7T ln 


\2 ) 

10 

LwiW 2 J J 


( 21 ) 


where a\ = 0.1381? 1 and l c is a constant, Q = w(l — cos[6b]/2)(l — cos[0 2 ]/2) is the modified frequency, 
and u>i = 2nu c Uj (1 — u c M 3 cos 0i) 1 ( 2/2 — 2/i) — 1 is the peak frequency. Here is the observer angle from 
local shock wave shear layer interactions to observer i. The first term of Eqn. 21 involving Dj enforces the 
growth of the length scale as the flow develops from the nozzle exit, and the second term contains frequency 
dependence. Frequency dependence of the length scale is inspired by the work of Lieb and Goldstein 34 and 
Morris and Boluriaan. 35 Separation between shocks within the second term is also included to emphasize 
the importance that large spatially coherent turbulence is responsible for BBSAN generation. The radial 
turbulent length scale is one-third the axial length scale. The time scale of turbulence, r s , is modeled as 


T s = 


1~cls 

U C U 


tanh 


4 u 

(Wi + U! 2 )_ 


( 22 ) 


where r c is a constant. Here, r s also contains frequency dependence and shock separation distance within 
its second factor and is consistent with Eqn. 21. The variation of turbulent kinetic energy, k, is modeled as 
a log-normal distribution 


k = a^u 2 exp 


-a 0 1 In 


Vi 

_Vc_ 


(23) 


and is calibrated empirically relative to the measurements of Bridges and Wernet, 10 where a 2 = 1 and 
a 3 = 11/300 are constants. Unfortunately, Eqn. 23 does not account for the large increases of k near each 
shock wave observed in measurements of off-design jets. 

Now that the arguments of the model involving an equivalent meanflow and basic turbulence statistics 
have been proposed, we can form a model for the coefficient matrix Aiji m 


A-ijim — Df exp 


— (^i — 0 2 ) 2 < 1 + St 


4/3 u c (l — u c Mj cos di) 


100 


Dj A 


x exp 


10 , 


— (6»i - 9 2 ) 2 St 1/2 - + - tanh[3(0 2 - tt/3)] 


x exp 


-[ls(u ~ u 2 ) + 4 iu(y 2 -yi)(u- &)] 
4u 2 


(24) 


where Pf is a constant, St is the Strouhal number based on u 3 and Dj, a tilde denotes that the modified 
frequency is used, and A is the axial distance between the first and second shock wave. Equation 24 is in 
a similar form to that proposed by Miller 26 for mixing noise cross-spectra, and its justification is discussed 
there. It contains two (relatively weak) terms that are dependent on the radiation directions and frequency. 

The prediction of spatial coherence is performed by evaluating 


T(xi,x 2 ,uj) 


G{x 1 ,x 2 ,u)G*{x 1 ,x 2 ,u;) 
G(xi, x 1 ,gj)G(x 2 , x 2 ,uj) ’ 


(25) 


where the superscript * represents the complex conjugate. The evaluation of T requires three separate 
evaluations of Eqn. 15. Equation 15 can now be evaluated using these simple models based on experimentally 
or numerically observed trends. 


Results 

This section shows comparisons of predicted auto-spectra and coherence of BBSAN using Eqn. 15 with 
various measurements. Comparisons of auto-spectra and spatial coherence are shown in the near-field and 
far-field. Measurements are obtained from various facilities, and their diversity demonstrates the robustness 
of the method. 
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Far-Field Auto-Spectra 

Predictions using Eqn. 15 for power spectral density (PSD) in the far-field are compared with measurement. 
Comparisons are shown on a non-dimensional basis using St and adjusting SPL to unit St via addition of 
101og 10 \ujDj 4 ] to the PSD. Measurements are obtained from Bridges and Brown 9 that were conducted at 
the Small Hot Jet Acoustic Rig (SHJAR) facility of NASA. We examine measurements from two axisymmetric 
nozzles from the SHJAR database. The first nozzle is conical and convergent with Aid of unity, and the 
second nozzle is designed with the method of characteristics with Aid = 1-50. Both nozzles have an exit 
diameter of D = 0.0508 m. Both over- and under-expanded, and heated and cold, conditions are shown. 
Microphones are located on a circular arc centered about the nozzle exit and measurements are corrected 
to a distance of R— lOOD by accounting for spherical spreading. Comparisons are conducted on a lossless 
basis, and measurements account for associated atmospheric absorption effects. 

Eight far-held conditions are considered. The first is shown in Fig. 1 for the convergent nozzle operating 
at Mj = 1.24 and TTR = 1.00. Solid lines represent the measurement from the SHJAR facility of Bridges 
and Brown, 9 and the dashed lines represent the predictions from Eqn. 15. The x-axis is St and the y-axis is 
SPL per unit St. There are five sets of comparisons at observer angles 4/ = 50 through 130 degs. separated 
by 20 deg. increments. Note that the y- axis lines are 30 AdB apart and the maximum broadband SPL per 
unit St level is written next to each measurement. The observer angle is measured from the upstream nozzle 
axis. Recall that the measurement represents the total noise and the prediction is only for the BBSAN 
component of the total noise spectrum. Predictions are in agreement with measurement with respect to 
the maximum BBSAN intensity and peak frequency. Multiple broad lobes are apparent in the upstream 
and sideline observer directions, and in the downstream direction broaden significantly. At frequencies lower 
than St ~ 0.3 the predictions fall-off quickly as the frequency bands are dominated by noise from turbulent 
mixing. At higher frequencies the predictions fall-off from measurement, and the reason for this is discussed 
below. 

We now examine the same moderately off-design condition from the convergent nozzle operating at 
Alj = 1.24 and increase the heating to TTR = 3.20. Predictions are compared with measurement in Fig. 2. 
Relative to Fig. 1, it is apparent that the addition of heating has eliminated the many shock-associated tones. 
In actuality the BBSAN intensity is relatively insensitive to changes in temperature as shown recently by 
Viswanathan et al., 3 Kuo et ah, 12 and Miller. 4 Predictions show the same level of agreement as the unheated 
case, where the peak intensity and peak frequency match measurement. 

Two highly off-design comparisons are shown in Figs. 3 and 4. The convergent nozzle operates at 
Alj = 1.47, and TTR = 1.00 or TTR = 3.20. Screech tones are present within both measurements. The 
unheated case in Fig. 3 shows relatively less agreement than the heated case. This is partly due to the effect 
of heating on the acoustic feedback loop within the jet plume. In all directions the peak predicted intensities 
are lower than measurement. This is not the case for the heated comparison at Alj = 1.47. 

A final heated comparison for the convergent nozzle is shown in Fig. 5. The jet operates at Alj = 1.57 
and TTR = 3.20. The spectrum at the sideline location is used to calibrate the coefficients present within 
Eqn. 15. Once the coefficients are calibrated, they are held constant irrespective of the jet operating condi- 
tions. It comes as no surprise that the spectral shape of the main BBSAN lobe matches the measurement. 
No optimization algorithm is used for this task, and certainly more optimized values can be chosen. In 
the upstream direction the peak intensity and peak frequency continue to agree with measurement. Un- 
fortunately, in the downstream radiation direction the spectral width of the BBSAN does not increase, as 
observed for select cases. It is emphasized that an optimized set of coefficients would increase the agreement 
for all predictions, but that is not the purpose of this investigation. 

We now turn our attention to the convergent-divergent method of characteristics nozzle with Aid = 1.50. 
Figure 6 shows comparisons of predictions using Eqn. 15 with measurements from the nozzle operating at 
Mj = 1.294 and TTR = 1.00. This flow is over-expanded unlike the others considered. The over-pressure of 
the tonal components is larger than other cases examined, and this highly affects the BBSAN component. 
The predicted low and high frequency BBSAN fall-off is too small and large, respectively. 

Two final far-held auto-spectral comparisons are shown in Figs. 7 and 8. The Md = 1.50 nozzle operates 
under-expanded at Alj = 1.70 and M ? = 1.828. Both hows are unheated. In general the overall trends are 
captured. It is interesting that for the Mj = 1.70 case the higher frequency components are captured by the 
model. Otherwise, the same trends are observed in these comparisons as those previously examined. 

The auto-spectral predictions shown here often exhibit a higher rate of intensity fall-off relative to mea- 
surement. Particular predictions of Morris and Miller 23 also showed this behavior to a lesser extent. The 
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method of Tam 20,22 used an empirical correction factor to control the intensity decay with frequency at high 
frequencies. Harper-Bourne and Fisher 1 avoided this problem by using a master source spectrum that was 
calibrated with measurement. An empirical correction, such as the one developed by Tam, 20:22 can easily 
be used to correct the high frequency intensity decay. 

The predicted peak intensities at the BBSAN peak frequency compared very favorably across the jet 
operating conditions examined at the sideline location. Refraction effects within the jet shear layer have 
little effect on the intensity of waves propagating to the far-held in the sideline direction. These observations 
give credibility to the equivalent source model of Eqn. 7. This choice of equivalent source combined with 
other model arguments allows for the majority of the characteristics of BBSAN to be predicted satisfactorily. 

Near-Field Auto-Spectra 

We now turn our attention to examining the characteristics of the BBSAN intensity in the jet near-held. 
For this task we use the excellent measurements of Yu. 5 Yu 5 conducted near-held measurements using a 
convergent-divergent nozzle with exit diameter of D = 0.0508 nr. The jet operates at Mj = 1.67 and TTR 
= 1.00. Coincidently, this operating condition is almost identical to that presented for far-held auto-spectral 
comparisons as shown in Fig. 7. Yu 5 used a near-held microphone array that transversed in the axial and 
radial directions over a very large spatial range. Measurements were processed by Yu in one-third octave 
center bands and presented as contour plots. Spatial dimensions are normalized by the nozzle exit diameter 
and an off-set in the radial direction of one-half D. It must be emphasized that the prediction model is for 
BBSAN and not shock-associated tones or turbulent mixing noise. It is expected that there will be some 
differences between prediction and measurement for this reason. 

Figure 9 shows near-held contours of SPL in the 16 kHz one-third octave center band from the Mj = 1.67 
off-design unheated jet. Contours are labeled by the local sound pressure level. The predicted levels are shown 
in Fig. 9(a) and the measured levels are shown in Fig. 9(b). The broad black line within the measurement 
was used by Yu to illustrate the region where microphones were not present. At this relatively low frequency 
(for a laboratory scale nozzle) a single broad lobe appears within both the predicted and measured contours 
that is directed slightly upstream relative to the sideline direction. Very close to the shock wave shear layer 
interactions (0 < x/D < 10 and 0 < r < D/2) the predictions show wavefronts surrounding the sources. 
Predictions have high enough resolution to capture the very near-field behavior while measurements lack 
detail. The agreement of overall levels and directivity at this frequency is fairly good, and as an example the 
contour level of 132 dB occurs both in prediction and measurement near x/D ss 2.5 and r/D — 1/2 ss 15. 

The same jet near-field is now examined at the 25 kHz one-third octave center band. These comparisons 
are shown in Fig. 10. Here, a single broad lobe is again present and its directivity remains close to the 
sideline direction. Relative to the previous comparison, the prediction is overall slightly less intense and 
points slightly more downstream. The measured contours of Yu 5 have a considerably thiner width, and the 
predicted contours are upstream by a few diameters. 

Figure 11 shows SPL contours at 31.5 kHz one-third octave center band. At this higher frequency we 
start to observe the formation of a second contour lobe in the upstream direction that is not present at 
lower frequencies. The more intense lobe has a directivity at a larger observer angle. Compared to the 
measurement the predicted intensity is slightly farther upstream and radiates at a slightly lower intensity. 

A final comparison of contours of SPL and measurements of Yu 5 is shown in Fig. 12. At the 40 kHz 
one-third octave band frequency the dominant radiation direction of both lobes is shifted to larger observer 
angles. The measurement also shows stronger double lobes. Contour levels between the two maps differ by 
a couple one-third octave dB and show some spatial shift. 

As an exercise we examine the region very close to the shock wave shear layer interactions of the off-design 
jet of Yu. 5 The M ( i = 1.50 nozzle operates at Mj = 1.67 and TTR = 1.00. Figure 13 shows contours of 
BBSAN intensity at one-third octave bands of 16, 25, 31.5, and 40 kHz. The contour maps are consistent 
between figures. Regions of high intensity near the nozzle lip line are very close to the shockwave shear layer 
interactions. By examining the contours near each shock wave shear layer interaction, it is clear that the 
complexity of the flow-field increases with increasing frequency. Finally, an important point of this theory 
is the highest intensity sources are located at multiple shocks downstream from the nozzle exit, and this 
is corroborated by other measurements that used beam-forming. This is mainly due to the growth of the 
intensity of turbulent kinetic energy with increasing axial distance from the nozzle exit and the turbulent 
kinetic energy interactions with the shock waves. 
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Far-Field Coherence 

Predicted spatial coherence in the jet far-held is now compared with measurements of Viswanathan et al. 3 
We focus on the comparison of coherence due to its relevance to practical problems. Predicted coherence 
is governed by Eqn. 25. The arguments of Eqn. 25 are provided by the developed model Eqn. 15. Note 
that to evaluate coherence between two observer positions (aq and aq) the cross-spectral prediction using 
Eqn. 25 is evaluated three times. Coherence is examined by holding the observer aq constant while varying 
the azimuthal position of observer X 2 . Reversing the two observer positions results in the same coherence. 

Figure 14 shows predicted coherence compared with measurements of Viswanathan et al. 3 These measure- 
ments were conducted by Viswanathan et al. 3 using a convergent nozzle with exit diameter of D = 0.06223 
m. Multiple microphones were located on an azimuthal arc in ten degree increments at the sideline direction. 
The nozzle operated at Mj = 1.57 and TTR = 3.20. The solid lines within Fig. 14 represent the coherence 
calculated by Viswanathan et al. 3 and the dashed lines are predictions using Eqn. 25. Four variations of 
azimuthal angle separation are shown at A cf> = 10, 20, 30, and 40 degs. As the azimuthal angle increases 
the coherence diminishes from a maximum of approximately 0.90 to 0.46. For this particular jet flow the 
high levels of coherence from turbulent mixing are present in the range of 0.001 ^ St 0.1. Coherence due 
to the BBSAN source is approximately in the range of 0.1 ;$ St ;$ 1.0. Coherence at higher frequencies is 
almost negligible for this particular jet operating condition in the far-held; however, for other jet operating 
conditions significant coherence might be present due to mixing noise or possibly shock-associated noise. 
The predicted maximum amplitude of coherence and peak frequency of coherence are in close agreement 
with measurements. The predicted envelope of coherence is perhaps too wide relative to measurement. 

The same nozzle and microphone array of Viswanathan et al. 3 are used to make a second comparison 
shown in Fig. 15. The convergent nozzle operates at Mj = 1.24 and TTR = 3.20. At this more conservative 
off-design condition the azimuthal coherence variation follows the same general trend as shown in Fig. 14. 
The spatially varying coherence due to the off-design condition is apparent as a broad lobe near the BBSAN 
peak frequency and intensity. Predicted peak coherence is within approximately 0.1 of measurement for the 
four angles examined. 

Near-Field Coherence 

Final comparisons are shown for spatial coherence in the jet near-field using Eqns. 15 and 25. Because the 
model is developed for the near-field, no alteration of the equations or their evaluation is required. Figure 16 
shows comparisons of predicted coherence with the measurements of Savarese et al. 13 Measurements were 
conducted with a convergent contoured nozzle with exit diameter of 0.04 nr operating at Mj = 1.22 and TTR 
= 1.00. A linear microphone array was located at R/D = 4 radially from the nozzle centerline axis and was 
parallel to the nozzle centerline axis. Each microphone was separated by a distance of 0.625.D. The reference 
microphone (aq) is located at the sideline location and R/D = 4. Observer positions of aq are located at 
Ax/D = -1.250, -0.625, 0.625, and 1.250. Predictions and measurements follow the same trends and show 
the same agreement as the far-held comparisons. In the near-held, there is much more variation in coherence, 
especially at higher frequencies, and a strong screech tone is present near St ss 0.45. At frequencies higher 
than St ~ 1, there is a combination of coherent nozzle tones (as discussed by Suzuki and Colonius 36 ) and 
some coherence due to BBSAN. The maximum predicted coherence matches the measured maximum near 
the BBSAN peak coherent frequency fairly well except for Ax/D = 0.625, where it is lower by approximately 
0.1. The peak frequencies are well predicted; that is to be expected relative to the far-held results. 

Two additional comparisons are shown using the predictions of Eqn. 15 and the measurements of 
Savarese et al. 13 These are presented in Figs. 17 and 18. Fully-expanded Mach numbers are 1.16 and 
1.28 respectively, and the other jet conditions and observer locations remain the same. Like the previous 
near-held case, as the axial distance increases the coherence magnitude decreases. The width of the BBSAN 
coherence generally captures the main BBSAN lobe. Higher and lower frequency coherence is likely due to 
other sources due to their dominance of the BBSAN within the pressure time history. 

Comparisons of predictions of acoustic intensity or coherence in the near-held are difficult given the 
large number of aerodynamic effects. The measurements are equally difficult. Because of the agreement of 
intensity in the far-held and near-held and coherence in the far-held, some conhdence is obtained for the 
prediction of coherence due to shock wave shear layer interactions in the near-held. The relative maximum 
BBSAN coherence and peak frequency in the near-held seems to be correctly predicted, though there is 
additional coherence at St > 1 that is likely due to aerodynamic effects or other noise sources. 
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Conclusion 


The cross-spectral acoustic analogy is used to predict auto-spectra and coherence of broadband shock- 
associated noise (based on cross-spectra) within the jet near-field and far-held. An equivalent source term is 
proposed for the shock-associated noise source that is consistent with other models for BBSAN. Arguments 
of the model involve flow-held statistics near the shock wave shear layer interactions. These arguments are 
dependent on the jet Mach number, temperature ratio, and nozzle exit diameter. A simplihed model is 
developed based on the assumption that the how is axisymmetric, and its arguments are dependent on a 
variety of measurements. The resultant model consists of a double summation of the source over all shock 
wave shear layer interactions. 

It is shown that predictions of auto-spectra and coherence generally compare favorably for a wide range of 
jet Mach numbers and temperature ratios. This observation gives credibility to the choice of the equivalent 
source model, which is satisfying given that refraction effects are not explicitly accounted for within this 
approach. Auto-spectral predictions exhibit a high rate of intensity fall-off relative to measurement, and 
this can easily be altered with an empirical correction factor. Near-held auto-spectral predictions correctly 
capture the formation of increasing numbers of directional lobes that originate from the shock cell structure. 
There are large variations of intensity between each shock wave shear layer interaction as observed in the am- 
plified view of the contours. Though, there is some axial shift within the predictions relative to measurement, 
their intensities are relatively similar. The complexity of coalescence of radiating broadband shock-associated 
noise increases greatly with increasing frequency. High intensity sources are located at multiple shock wave 
shear layer interactions downstream from the nozzle exit, and this is corroborated with measurements using 
beam-forming. The log-normal model for turbulent kinetic energy is essential for capturing this trend. 

At mid-frequencies within the near-field or far-field, strong broad lobes of coherence are due to shock 
wave shear layer interactions, while at lower and higher frequencies this coherence is due to other effects. 
Measured coherence in the near-field contains additional near-field effects due to other sources and additional 
hydrodynamic fluctuations. Unfortunately, in some cases the predicted envelope of the coherence is too wide 
relative to measurement. Predicted near-field coherence is nearly zero at high frequencies, and measurement 
shows coherent tones dominating high frequencies. It is likely that the spatial BBSAN coherence is nearly 
zero at high frequencies based upon the far-field comparisons, where strongly coherent nozzle tones are 
not present. Given the wide range of nozzle operating conditions and the difficulty of the measurements, 
predictions capture the overall trends of the statistics of the near-field shock-associated noise satisfactorily. 
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Figure 1. Predictions of BBS AN intensity compared 

with measurement for Md — 1.00, Mj — 1.24, 
and TTR = 1.00. 



Figure 3. Predictions of BBSAN intensity compared 
with measurement for Md — 1.00, Mj = 1.47, 
and TTR = 1.00. 



Figure 2. Predictions of BBSAN intensity compared 
with measurement for Md — 1.00, Mj — 1.24, 
and TTR = 3.20. 



Figure 4. Predictions of BBSAN intensity compared 
with measurement for Md — 1.00, Mj — 1.47, 
and TTR = 3.20. 

14 of 21 


American Institute of Aeronautics and Astronautics 



Figure 5. Predictions of BBSAN intensity compared 
with measurement for Md — 1 . 00 , Mj — 1 . 57 , 
and TTR = 3.20. 



Figure 7. Predictions of BBSAN intensity compared 
with measurement for Md — 1 . 50 , Mj = 1 . 70 , 
and TTR = 1.00. 



Figure 6. Predictions of BBSAN intensity compared 
with measurement for Md — 1 . 50 , Mj — 1 . 294 , 
and TTR = 1.00. 



Figure 8. Predictions of BBSAN intensity compared 
with measurement for Md — 1 . 50 , Mj — 1 . 828 , 
and TTR = 1.00. 
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(a) BBSAN Prediction (b) Measurement 

Figure 9. Near-field contours of SPL at 16 kHz. The M d = 1.50 nozzle operates at Mj = 1.67 and TTFt = 1.00. 




(a) BBSAN Prediction 


(b) Measurement 


Figure 10. Near-field contours of SPL at 25 kHz. The = 1.50 nozzle operates at Mj = 1.67 and TTR = 1.00. 
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(a) BBSAN Prediction 

Figure 11. Near-field contours of SPL at 31.5 kHz. The 



x/D 

(b) Measurement 

= 1.50 nozzle operates at Mj = 1.67 and TTR = 1.00. 




Figure 12. 


(a) BBSAN prediction 

Near-field contours of SPL at 40.0 kHz. 


(b) Measurement 

The Md = 1.50 nozzle operates at Mj = 1.67 and TTR. = 1.00. 
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r/D -1/2 r/D- 1/2 




(a) 16 kHz (b) 25 kHz 




(c) 31.5 kHz (d) 40 kHz 

Figure 13. Predicted near-field contours of BBSAN intensity. The Md = 1.50 nozzle operates at Mj = 1.67 and TTR 

= 1 . 00 . 
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Figure 14. Predicted coherence compared with measurement for Md = 1 . 00 , Mj = 1 . 57 , and TTR = 3.20. Comparisons 
are in the far-field at the sideline and R/D = 100 . The observer microphone varies azimuthally. 



Figure 15. Predicted coherence compared with measurement for Md = 100 , Mj = 1 . 24 , and TTR = 3.20. Comparisons 
are in the far-field at the sideline and R/D — 100 . The observer microphone varies azimuthally. 
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Figure 16. Predicted coherence compared with measurement for Md = 1.00, Mj = 1.22, and TTR = 1 . 00 . Comparisons 
are in the near-field at the sideline and R/D = 4. The observer microphone varies axially. 






Figure IT. Predicted coherence compared with measurement for Md = 1.00, Mj = 1.16, and TTR. = 1 . 00 . Comparisons 
are in the near-field at the sideline and R/D = 4. The observer microphone varies axially. 
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Figure 18 . Predicted coherence compared with measurement for Md = 1.00, Mj = 1.28, and TTR = 1 . 00 . Comparisons 
are in the near-field at the sideline and R/D = 4. The observer microphone varies axially. 
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